Carbon stock projection for four major forest plantation species in Japan.

Carbon sequestration via afforestation and forest growth is effective for mitigating global warming. Accurate and robust information on forest growth characteristics by tree species, region, and large-scale land-use change is vital and future prediction of forest carbon stocks based on this information is of great significance. These predictions allow exploring forestry practices that maximize carbon sequestration by forests, including wood production. Forest inventories based on field measurements are considered the most accurate method for estimating forest carbon stocks. Japan's national forest inventories (NFIs) provide stand volumes for all Japanese forests, and estimates from direct field observations (m-NFIs) are the most reliable. Therefore, using the m-NFI from 2009 to 2013, we selected four major forest plantation species in Japan: Cryptomeria japonica, Chamaecyparis obtusa, Pinus spp., and Larix kaempferi and presented their forest age-carbon density function. We then estimated changes in forest carbon stocks from the past to the present using the functions. Next, we investigated the differences in the carbon sequestration potential of forests, including wood production, between five forestry practice scenarios with varying harvesting and afforestation rates, until 2061. Our results indicate that, for all four forest types, the estimates of growth rates and past forest carbon stocks in this study were higher than those considered until now. The predicted carbon sequestration from 2011 to 2061, assuming that 100 % of harvested carbon is retained for a long time, twice the rate of harvesting compared to the current rate, and a 100 % afforestation rate in harvested area, was three to four times higher than that in a scenario with no harvesting or replanting. Our results suggest that planted Japanese forests can exhibit a high carbon sequestration potential under the premise of active management, harvesting, afforestation, and prolonging the residence time of stored carbon in wood products with technology development.

Cambium Reactivation Is Closely Related to the Cell-Cycle Gene Configuration in Larix kaempferi.

Dormancy release and reactivation in temperate trees are mainly controlled by temperature and are affected by age, but the underlying molecular mechanisms are still unclear. In this study, we explored the effects of low temperatures in winter and warm temperatures in spring on dormancy release and reactivation in Larix kaempferi. Further, we established the relationships between cell-cycle genes and cambium cell division. The results showed that chilling accelerated L. kaempferi bud break overall, and the longer the duration of chilling is, the shorter the bud break time is. After dormancy release, warm temperatures induced cell-cycle gene expression; when the configuration value of the cell-cycle genes reached 4.97, the cambium cells divided and L. kaempferi reactivated. This study helps to predict the impact of climate change on wood production and provides technical support for seedling cultivation in greenhouses.

Stem moisture content prediction model for Larix olgensis based on beta regression.

To investigate the longitudinal variation patterns of sapwood, heartwood, bark and stem moisture content along the trunk of artificial Larix olgensis, we constructed mixed effect models of moisture content based on beta regression by combining the effects of sampling plot and sample trees. We used two sampling schemes to calibrate the model, without limiting the relative height (Scheme Ⅰ) and with a limiting height of less than 2 m (Scheme II). The results showed that sapwood and stem moisture content increased gradually along the trunk, heartwood moisture content decreased slightly and then increased along the trunk, and bark moisture content increased along the trunk and then levelled off before increasing. Relative height, height to crown base, stand area at breast height per hectare, age, and stand dominant height were main factors driving moisture content of L. olgensis. Scheme Ⅰ showed the stable prediction accuracy when randomly sampling moisture content measurements from 2-3 discs to calibrate the model, with the mean absolute percentage error (MAPE) of up to 7.2% for stem moisture content (randomly selected 2 discs), and the MAPE of up to 7.4%, 10.5% and 10.5% for sapwood, heartwood and bark moisture content (randomly selected 3 discs), respectively. Scheme Ⅱ was appropriate when sampling moisture content measurements from discs of 1.3 and 2 m height and the MAPE of sapwood, heartwood, bark and stem moisture content reached 7.8%, 11.0%, 10.4% and 7.1%, respectively. The prediction accuracies of all mixed effect beta regression models were better than the base model. The two-level mixed effect beta regression models, considering both plot effect and tree effect, would be suitable for predicting moisture content of each part of L. olgensis well.

Phylogeny and expression patterns of ERF genes that are potential reproductive inducers in hybrid larch.

BACKGROUND: Larch is an important component of northern forests and a major cultivated tree species in restoration of forest cover using improved seed material. In recent years, the continuous low seed production has severely affected the production of improved variety seedlings and natural regeneration. However, research on the reproductive growth of gymnosperms is extremely scarce. RESULTS: In this study, based on differential transcriptome analysis of two asexual reproductive phases, namely high-yield and low-yield, we further screened 5 ERF family genes that may affect the reproductive development of larch. We analyzed their genetic relationships and predicted their physicochemical properties. The expression patterns of these genes were analyzed in different tissues, developmental stages, hormone treatments, and environmental conditions in hybrid larch. CONCLUSION: The results showed that all 5 genes were induced by low temperature and ABA, and their expression patterns in different tissues suggested a suppressive role in the development of female cones in larch. Among them, LkoERF3-like1 and LkoERF071 may be involved in the flowering age pathway. This study enriches the scarce research on reproductive development in gymnosperms and provides a theoretical basis and research direction for regulating the reproductive development of larch in seed orchards.

A new perspective on structural characterisation and immunomodulatory activity of arabinogalactan in Larix kaempferi from Qinling Mountains.

In this study, crude polysaccharide (LAG-C) and homogeneous arabinogalactan (LAG-W) were isolated from Qinling Larix kaempferi of Shaanxi Province. Bioactivity assays showed that LAG-W and LAG-C enhanced the phagocytic ability, NO secretion, acid phosphatase activity, and cytokine production (IL-6, IL-1ß, and TNF-α) of RAW264.7 macrophages. Notably, LAG-W exhibited a significantly stronger immunomodulatory effect than LAG-C. The primary structure of LAG-W was characterised by chemical methods (monosaccharide composition, methylation analysis, and alkali treatment) and spectroscopic techniques (gas chromatography-mass spectrometry, high-performance liquid chromatography-mass spectrometry, and 1D/2D nuclear magnetic resonance). LAG-W was identified as a 22.08 kilodaltons (kDa) neutral polysaccharide composed of arabinose and galactose at a 1:7.5 molar ratio. Its backbone consisted of repeated →3)-ß-Galp-(1→ residues. Side chains, connected at the O-6 position, were mainly composed of T-ß-Galp-(1→ and T-ß-Galp-(1→6)-ß-Galp-(1→ residues. And it also contained small amounts of T-ß-Arap-(1→, T-α-Araf-(1→6)-ß-Galp-(1→6)-ß-Galp-(1→, and T-α-Araf-(1→3)-α-Araf-(1→6)-ß-Galp-(1→ residues. By structurally and functionally characterising L. kaempferi polysaccharides, this study opens the way for the valorisation of this species.

Photosynthetic responses of Larix kaempferi and Pinus densiflora seedlings are affected by summer extreme heat rather than by extreme precipitation.

The frequency and intensity of summer extreme climate events are increasing over time, and have a substantial negative effect on plants, which may be evident in their impact on photosynthesis. Here, we examined the photosynthetic responses of Larix kaempferi and Pinus densiflora seedlings to extreme heat (+ 3 °C and + 6 °C), drought, and heavy rainfall by conducting an open-field multifactor experiment. Leaf gas exchange in L. kaempferi showed a decreasing trend under increasing temperature, showing a reduction in the stomatal conductance, transpiration rate, and net photosynthetic rate by 135.2%, 102.3%, and 24.8%, respectively, in the + 6 °C treatment compared to those in the control. In contrast, P. densiflora exhibited a peak function in the stomatal conductance and transpiration rate under + 3 °C treatment. Furthermore, both species exhibited increased total chlorophyll contents under extreme heat conditions. However, extreme precipitation had no marked effect on photosynthetic activities, given the overall favorable water availability for plants. These results indicate that while extreme heat generally reduces photosynthesis by triggering stomatal closure under high vapor pressure deficit, plants employ diverse stomatal strategies in response to increasing temperature, which vary among species. Our findings contribute to the understanding of mechanisms underlying the photosynthetic responses of conifer seedlings to summer extreme climate events.

Distribution of exchangeable Ca2+ during the process of Larix decidua Mill. pollination and germination.

The involvement of Ca2+ ions in angiosperms sexual processes is well established, while in gymnosperms, such knowledge remains limited and is still a topic of discussion. In this study, we focused on Larix decidua, using Alizarin-red S staining and the pyroantimonate method to examine the tissue and subcellular distribution of free and loosely bound Ca2+ ions at different stages of the male gametophyte's development and its interaction with the ovule. Our findings show that in larch, both the germination of pollen grains and the growth of pollen tubes occur in an environment rich in Ca2+. These ions play a crucial role in the adhesion of the pollen grain to the stigmatic tip and its subsequent movement to the micropylar canal. There is a significant presence of free and loosely bound Ca2+ ions in both the fluid of the micropylar canal and the extracellular matrix of the nucellus. As the pollen tube extends through the nucellus, we observed a notable accumulation of Ca2+ ions just above the entry to the mature archegonium, a region likely crucial for the male gametophyte's directional growth. Meanwhile, the localized presence of free and loosely bound Ca2+ ions within the egg cell cytoplasm may inhibit the pollen tubes growth and rupture, playing an important role in fertilization.

Host-specific growth responses of Larix kaempferi and Quercus acutissima to Asian gypsy moth defoliation in central Korea.

As the risk of gypsy moth outbreaks that have detrimental effects on forest ecosystem in the Northern Hemisphere increase due to climate change, a quantitative evaluation of the impact of gypsy moth defoliation is needed to support the adaptive forest management. To evaluate the host-specific impact of gypsy moth defoliation, radial growth and annual carbon accumulation were compared for one severely defoliated (Larix kaempferi (Lamb.) Carrière) and one moderate defoliated (Quercus acutissima Carruth.) host, in defoliated and non-defoliated site using tree-ring analysis. Finally, the resilience indices of radial growth variables were calculated to assess the ability of sampled trees to withstand defoliation. Gypsy moth defoliation mainly decreased latewood width and caused reduction in annual carbon absorption more than 40% for both tree species. However, L. kaempferi, showed the reduced growth until the year following defoliation, while Q. acutissima, showed no lagged growth depression and rapid growth recover. The findings show how each species reacts differently to gypsy moth defoliation and highlight the need of managing forests in a way that takes resilient tree species into account.

H2O2 Significantly Affects Larix kaempferi × Larix olgensis Somatic Embryogenesis.

Larch is widely distributed throughout the world and is an important species for timber supply and the extraction of industrial raw materials. In recent years, the hybrid breeding of Larix kaempferi and Larix olgensis has shown obvious heterosis in quick-growth, stress resistance and wood properties. However, its growth and development cycle is too long to meet general production needs. In order to shorten the breeding cycle, we have for the first time successfully established and optimized a somatic embryogenesis system for Larix kaempferi × Larix olgensis. We found that the highest rate of embryonal-suspensor mass (ESM) induction was observed when late cotyledonary embryos were used as explants. The induced ESMs were subjected to stable proliferation, after which abscisic acid (ABA) and polyethylene glycol (PEG) were added to successfully induce somatic embryos. Treatment with PEG and ABA was of great importance to somatic embryo formation and complemented each other's effect. ABA assisted embryo growth, whereas PEG facilitated the formation of proembryo-like structures. On top of this, we studied in more detail the relationship between redox homeostasis and the efficiency of somatic embryogenesis (frequency of ESM induction). During subculture, we observed the gradual formation of three distinct types of ESM. The Type I ESM is readily able to form somatic embryos. In contrast to type I, the type III ESM suffers from severe browning, contains a higher level of hydrogen peroxide (H2O2) and demonstrates a decreased ability to form somatic embryos. External treatment with H2O2 decreased the somatic embryogenesis efficiency of Type I and type III ESMs, or the higher the exogenous H2O2 content, the lower the resulting somatic embryogenesis efficiency. We found that treatment with the H2O2 scavenger DMTU (dimethylthiourea) could significantly increase the somatic embryogenesis efficiency of the type III ESM, as a result of a decline in endogenous H2O2 content. Overall, these findings have contributed to setting up a successful somatic embryogenesis system for larch production.

Genomic dissection of additive and non-additive genetic effects and genomic prediction in an open-pollinated family test of Japanese larch.

Genomic dissection of genetic effects on desirable traits and the subsequent use of genomic selection hold great promise for accelerating the rate of genetic improvement of forest tree species. In this study, a total of 661 offspring trees from 66 open-pollinated families of Japanese larch (Larix kaempferi (Lam.) Carrière) were sampled at a test site. The contributions of additive and non-additive effects (dominance, imprinting and epistasis) were evaluated for nine valuable traits related to growth, wood physical and chemical properties, and competitive ability using three pedigree-based and four Genomics-based Best Linear Unbiased Predictions (GBLUP) models and used to determine the genetic model. The predictive ability (PA) of two genomic prediction methods, GBLUP and Reproducing Kernel Hilbert Spaces (RKHS), was compared. The traits could be classified into two types based on different quantitative genetic architectures: for type I, including wood chemical properties and Pilodyn penetration, additive effect is the main source of variation (38.20-67.46%); for type II, including growth, competitive ability and acoustic velocity, epistasis plays a significant role (50.76-91.26%). Dominance and imprinting showed low to moderate contributions (< 36.26%). GBLUP was more suitable for traits of type I (PAs = 0.37-0.39 vs. 0.14-0.25), and RKHS was more suitable for traits of type II (PAs = 0.23-0.37 vs. 0.07-0.23). Non-additive effects make no meaningful contribution to the enhancement of PA of GBLUP method for all traits. These findings enhance our current understanding of the architecture of quantitative traits and lay the foundation for the development of genomic selection strategies in Japanese larch.

Anatomical and morphological changes in Pinus sylvestris and Larix sibirica needles under impact of emissions from a large aluminum enterprise.

Species-specific anatomical and morphological characteristics of Pinus sylvestris and Larix sibirica needles were studied at different levels of tree stand pollution by aluminum smelter emissions. The anatomical characteristics of the needle were studied using light microscopy. The level of tree stand pollution was determined using the cluster analysis outcomes of the pollutant elements content (fluorine, sulfur, and heavy metals) in the needles. Four levels of tree stand pollution were separated: low, moderate, high, and critical, as well as background tree stand in unpolluted areas. It was found that the state of tree phytomass deteriorated with increasing levels of pollution (from low to critical): pine crown defoliation increased to 85%, and larch defoliation increased to 65%. The life span of pine needles was reduced to 2-3 years, with a background value of 6-7 years. The change of morphological parameters was more pronounced in P. sylvestris: the weight and length of the 2-year-old shoot decreased by 2.7-3.1 times compared to the background values; the weight of needles on the shoot and the number of needle pairs on the shoot-by 1.9-2.1 times. The length of the needle and shoot and the number of L. sibirica brachyblasts decreased by 1.8-1.9 times. The anatomical parameters of the needle also changed to a greater extent in P. sylvestris. Up to the high level of tree pollution, we observed a decrease in the cross-sectional area of the needle, central cylinder, vascular bundle, area and thickness of mesophyll, number and diameter of resin ducts by 18-66% compared to background values. At the critical pollution level, when the content of pollutant elements in pine needles reached maximum values, the anatomical parameters of the remaining few green needles were close to background values. In our opinion, this may be due to the activation of mechanisms aimed at maintaining the viability of trees. A reduction in thickness and area of assimilation tissue in the L. sibirica needle was detected only at the critical pollution level. An upward trend in these parameters was found at low, medium, and high pollution levels of tree stand, which may indicate an adaptive nature. The results suggested that at a similar pollution level of trees, the greatest amount of negative anatomical and morphological changes were recorded in pine needles, which indicates a greater sensitivity of this species to technogenic emissions.

Rapid migration of Mongolian oak into the southern Asian boreal forest.

The migration of trees induced by climatic warming has been observed at many alpine treelines and boreal-tundra ecotones, but the migration of temperate trees into southern boreal forest remains less well documented. We conducted a field investigation across an ecotone of temperate and boreal forests in northern Greater Khingan Mountains of northeast China. Our analysis demonstrates that Mongolian oak (Quercus mongolica), an important temperate tree species, has migrated rapidly into southern boreal forest in synchrony with significant climatic warming over the past century. The average rate of migration is estimated to be 12.0 ± 1.0 km decade-1 , being slightly slower than the movement of isotherms (14.7 ± 6.4 km decade-1 ). The migration rate of Mongolian oak is the highest observed among migratory temperate trees (average rate 4.0 ± 1.0 km decade-1 ) and significantly higher than the rates of tree migration at boreal-tundra ecotones (0.9 ± 0.4 km decade-1 ) and alpine treelines (0.004 ± 0.003 km decade-1 ). Compared with the coexisting dominant boreal tree species, Dahurian larch (Larix gmelinii), temperate Mongolian oak is observed to have significantly lower capacity for light acquisition, comparable water-use efficiency but stronger capacity to utilize nutrients especially the most limiting nutrient, nitrogen. In the context of climatic warming, and in addition to a high seed dispersal capacity and potential thermal niche differences, the advantage of nutrient utilization, reflected by foliar elementomes and stable nitrogen isotope ratios, is also likely a key mechanism for Mongolian oak to coexist with Dahurian larch and facilitate its migration toward boreal forest. These findings highlight a rapid deborealization of southern Asian boreal forest in response to climatic warming.

Higher latewood to earlywood ratio increases resistance of radial growth to severe droughts in larch.

As drought has caused great losses of tree growth across the world, the mechanism of how trees adapt to drought has been extensively investigated. However, how trees change their late- to earlywood ratio (LER) to adapt to severe drought events remains poorly understood. We used a network of Larix principis-rupprechtii earlywood and latewood width data from 1979 to 2018, covering most of the distribution of planted larch across North China, to investigate how latewood proportion affected trees' resistance to drought. The interactions among LER, minimum temperature, vapor pressure deficit (VPD), growing season length, and their contributions to drought resistant (Rt) were estimated using structural equation models. The results show a significant increase in LER of the juvenile wood throughout the first 15 growth rings after which it stabilizes. The LER decreased significantly with elevation for the juvenile wood. March-May temperature and VPD were the main determinant in the LER of mature wood. The sensitivity of radial growth to droughts was positively changed with LER when LER was below 0.50, but negatively changed with LER when LER is above 0.50. We confirmed that high LER increases resistance of tree growth to severe droughts in L. principis-rupprechtii. Our results highlight that a higher proportion of latewood is formed in dry years, and this high drought sensitivity of LER in turn led to an increased resistance to drought. This combination of reduced radial growth during dry years, while the latewood proportion remains increases maybe an adaptive strategy of larch trees to cope with severe droughts.

UAV Photogrammetry for Estimating Stand Parameters of an Old Japanese Larch Plantation Using Different Filtering Methods at Two Flight Altitudes.

Old plantations are iconic sites, and estimating stand parameters is crucial for valuation and management. This study aimed to estimate stand parameters of a 115-year-old Japanese larch (Larix kaempferi (Lamb.) Carrière) plantation at the University of Tokyo Hokkaido Forest (UTHF) in central Hokkaido, northern Japan, using unmanned aerial vehicle (UAV) photogrammetry. High-resolution RGB imagery was collected using a DJI Matrice 300 real-time kinematic (RTK) at altitudes of 80 and 120 m. Structure from motion (SfM) technology was applied to generate 3D point clouds and orthomosaics. We used different filtering methods, search radii, and window sizes for individual tree detection (ITD), and tree height (TH) and crown area (CA) were estimated from a canopy height model (CHM). Additionally, a freely available shiny R package (SRP) and manually digitalized CA were used. A multiple linear regression (MLR) model was used to estimate the diameter at breast height (DBH), stem volume (V), and carbon stock (CST). Higher accuracy was obtained for ITD (F-score: 0.8-0.87) and TH (R2: 0.76-0.77; RMSE: 1.45-1.55 m) than for other stand parameters. Overall, the flying altitude of the UAV and selected filtering methods influenced the success of stand parameter estimation in old-aged plantations, with the UAV at 80 m generating more accurate results for ITD, CA, and DBH, while the UAV at 120 m produced higher accuracy for TH, V, and CST with Gaussian and mean filtering.

Adaptation of Betula pendula Roth., Pinus sylvestris L., and Larix decidua Mill. to environmental stress caused by tailings waste highly contaminated by trace elements.

The seedlings of some tree species can successfully develop in areas polluted by heavy metals. Research on such species is important in order to explore the possibility of introducing tree species for the permanent biological stabilization and reclamation of post-flotation tailings, especially after the final recycling of trace metals, but where concentrations remain much higher than in natural soils. To better understand the adaptation and reaction of Betula pendula Roth., Pinus sylvestris L., and Larix decidua Mill. seedlings to heavy metals pollution caused by tailings waste highly contaminated by trace elements: 1) the relationships between the concentration of heavy metals in the soil substrate, the efficiency of heavy metal ions accumulation in plant organs, and the biometric parameters of the seedlings; and 2) the threshold content of heavy metals in the roots above which the plant physiological response is triggered was determined. We assume that there are certain limit concentrations of heavy metals in the soil and fine roots, which depend on the tree species and beyond which the plant responds strongly to stressThe obtained results showed that Betula is a suitable species for the phytostabilization of post-flotation tailings due to its rapid growth rate and production of root biomass. The accumulation of metals in Betula roots was found to be much greater than in Pinus and Larix. Despite the high concentrations of heavy metals in the prepared substrates, there was only a slight transfer of these elements to the aboveground parts of the plant. At high soil concentrations, the heavy metals adversely affected the cellular and physiological processes of plants. In plants growing in such conditions, the activity of the antioxidant system depended both on the species and organ of the plant, as well as on the type and metal concentration.

Construction of universal equations for knot attributes of three coniferous species.

We constructed base model, dummy variable model, and mixture model with three variables including knot diameter, loose knot length, and sound knot length with three typical coniferous species, Pinus koraiensis, Larix olgensis, and Pinus sylvestris var. mongolica, from the Linkou Forestry Bureau and Mengjiagang forest farm in Heilongjiang Province in 2020. We analyzed the differences in knot properties among different tree species and simplified the modeling work. Firstly, we collected relevant knot property data through the sectioning method based on relevant literature, transformation of the model form and substitution of related variables to conduct a base model. We transformed the species into dummy variables as qualitative factors, and introduced the dummy variable model of the relevant attributes into the base model. We introduced the random effects of sample trees and sample plots when constructing the mixture model. By comparing evaluation indicators, such as Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC), the mixture model with the best fitting effect was selected. We selected the optimal universal equation by comparing the fitting accuracy of the base model, dummy variable model and mixture model. The fitting accuracy of the dummy variable model and mixture model was higher than that of the basic model. The evaluation indicators (AIC and BIC) showed that the mixture model had a better fitting effect on knot properties than the dummy variable model. In the model comparison results, R2 of mixture models for sound knot length, the loose knot length, and knot diameter increased by 13.2%, 84.8% and 40.3%, respectively. The predictive accuracy of the three base models for different tree species' knot attributes was above 90%, and both the prediction accuracy of the dummy variable model and mixture model were above 94%, indicating that the constructed models could well predict knot-related properties. From the perspective of tree species, the sound knot length, knot diameter, and loose knot length was in order of P. sylvestris var. mongolica > P. koraiensis > L. olgensis. Fitted results of the dummy variable model and the mixture model were superior to the basic model, with higher accuracy.

Diversity and negative effect of PM0.3-10.0 adsorbed by needles of urban trees in Irkutsk, Russia.

The study was performed in natural forests preserved within the Boreal zone city, Irkutsk, Russia. Test sites were selected in the forests in different districts of the city, where samples of Scots pine (Pinus sylvestris L.) and Siberian larch (Larix sibirica Ledeb.) needles were taken to study the adsorption on their surface of aerosol particles of different sizes, in microns: PM0.3, PM0.5, PM1, PM2.5, PM5, PM10. Scanning electron microscopy was used to obtain high-resolution photographs (magnification 800- × 2000, × 16,000) and aerosol particles (particulate matter-PM) were shown to be intensively adsorbed by the surface of needles, with both size and shape of the particles characterized by a wide variety. Pine needles can be covered with particles of solid aerosol by 50-75%, stomata are often completely blocked. Larch needles often show areas, which are completely covered with aerosol particles, there are often found stomata deformed by the penetration of PMx. X-ray spectral microanalysis showed differences in the chemical composition of adsorbed PMx, the particles can be metallic if metals predominate in their composition, carbonaceous-in case of carbon predominance-or polyelemental if the composition is complex and includes significant quantities of other elements besides metals and carbon (calcium, magnesium, potassium, sodium, sulfur, chlorine, fluorine). Since the particles contain a large proportion of technogenic pollutants, accumulation by the needles of some widespread pollutants was investigated. A direct correlation of a highly significant level between the concentration of PMx in the air and the accumulation of many heavy metals in pine and larch needles, as well as sulfur, fluorine, and chlorine, has been revealed, which indicates a high cleaning capacity of urban forests. At the same time, the negative impact of PMx particles on the vital status of trees is great, which shows in intense disturbance of the parameters of photosynthesis and transpiration, leading to a significant decrease in the growth characteristics of trees and reduction in the photosynthetic volume of the crowns. We consider that the results obtained are instrumental in developing an approach to improvement of urban forests status and creating a comfortable urban environment for the population.

Compatibility predictive model for regeneration quantities of Larix gmelinii natural forest in Daxing'anling Mountains, China.

The natural regeneration grade is an important foundation for formulating forest management measures. Traditional studies have only considered the regeneration quantities predictive model of the total stand or dominant tree species, but the consistency among the prediction results of different tree species and the total regeneration quantities of stand is not solved. That is, the regeneration prediction results at the stand level are not equal to the sum of the predicted results of all tree species. To address this, on the basis of the traditional counting model, we attempted to construct a compatibility predictive model for regeneration quantities of different tree species within the stand, which would provide a theoretical basis for the rational management and decision-making of natural forest. Based on the survey data from 96 standard plots of Cuigang Forest Farm, Xinlin Forest Farm, and Zhuangzhi Forest Farm in Daxing'an Mountains, we selected 30 basic indices from five aspects of site factor, soil factor, stand factor, tree diversity and stand spatial structure, and used Poisson model and negative binomial model as the basic models to construct the regeneration prediction models of Larix gmelinii, Betula platyphylla and other tree species. By comparing the accuracy and fitting effect of the two traditional counting models, we selected the optimal model and used the seemingly unrelated regressions to further construct the compatibility predictive model for regeneration quantities of different tree species. Poisson model was the best one for the regeneration of L. gmelinii, B. platyphylla, and other tree species. The test index RMSE of the compatibility predictive model for regeneration quantities of L. gmelinii, B. platyphylla, other tree species and total stand regeneration quantities were 388, 413, 504, and 871 trees·hm-2, respectively. The adjusted R2 was 0.389, 0.421, 0.488, and 0.407, respectively. The most influential variables for regeneration quantities of L. gmelinii, B. platyphylla and other tree species were Pielou evenness index of DBH (25.2%), herbal coverage (34.6%) and organic matter in B layer (23.2%). In this study, the compatibility predictive model system for regeneration quantities satisfied the additive logic among L. gmelinii, B. platyphylla, other tree species, and total stands, and provided a basis for accurately estimating natural regeneration.

Height-diameter models of regenerated saplings of Larix gmelinii based on dummy variable and quantile regression.

Based on data collected from 2054 saplings of Larix gmelinii forest in 55 fixed plots in 2018-2019 in Cuigang Forestry Station, Daxing'anling area, we classified the stand density index (SDI) into four classes, i.e., Class Ⅰ (SDI1<1863 plants·hm-2), Class Ⅱ (1863 plants·hm-2≤SDI2<2155 plants·hm-2), Class Ⅲ (2155 plants·hm-2≤SDI3<2459 plants·hm-2) and Class Ⅳ (SDI4≥2459 plants·hm-2) by using the quartile method. We constructed a dummy variable model and quantile regression model for the height-breast diameter of saplings of L. gmelinii with dummy variable method introduced SDI. The results showed that among the five selected representative non-linear tree height curve models, the Richards model fitted the best, with Ra2, RMSE and MAE of 0.7637, 0.8250 m and 0.5696 m. The dummy variable model including the SDI constructed based on the Richards model showed a 1.3% increase in Ra2 compared with the base model, while RMSE, MAE, and AIC decreased by 2.1%, 1.5%, and 11.2%, respectively. When the quantile τ was 0.5, Ra2 of quantile regression model was the maximum, and RMSE, MAE, AIC was the minimum, being 0.7612, 0.8294 m, 0.5657 m, and -767.19, respectively. Compared with SDI1, sapling height in SDI2-SDI4 was increased by 5.6%, 5.6%, and 11.3%, suggesting reasonable that regulation of stand density was conducive to increase the height growth of saplings in regeneration.

[Responses of radial growth of Larix principis-rupprechtii to abrupt warming]. / 华北落叶松径向生长对升温突变的响应.

To investigate the impacts of abrupt warming on tree growth, we collected tree ring cores from larch (Larix principis-rupprechtii) in three different sites, including Saihanba National Nature Reserve, Pangquangou National Nature Reserve, and Fengning Qiansongba National Forest Park. Based on the tree ring method, Mann-Kendall test was used to examine the occurrence time of temperature rise mutation. We further analyzed the radial growth law of larch before and after the temperature mutation and its correlation with the monthly climate data. The results showed that the sudden temperature rise occurred in the Saihanba area in 1987, the Fengning area in 1989, and the Pangquangou area 1994. Before the sudden warming, there was no significant trend for the radial growth in all the three regions. After the sudden warming, however, it decreased significantly (with a decrease rate of 0.08·10 a-1) in Saihanba area. The radial growth of larch increased significantly in Pangquangou area (with an increase rate of 0.10·10 a-1), while no significant change was observed in the Fengning area. Before the sudden warming, there was a significant positive correlation between the radial growth of larch in the Saihanba area and the highest temperature in May and June. After the sudden warming, there was a significant positive correlation with precipita-tion in July, and a highly significant positive correlation with the standardized precipitation evapotranspiration index(SPEI) from September of the previous year to July. Prior to the sudden warming, there was no significant relationship between the radial growth of larch in the Pangquangou area and monthly climate factors. However, after the sudden warming, a significant positive correlation was found with the lowest temperature in September of the pre-vious year. Before the sudden warming, the radial growth of larch in Fengning area was significantly negatively correlated with the lowest average temperature in July. After the sudden warming, it showed a significant negative correlation with the average and highest temperatures in June. Accordingly, the radial growth of larch in the Saihanba area experienced drought stress following a sudden temperature change. If temperature continues to rise in the future, larch in the Fengning area would also face drought stress. Conversely, warming conditions would be beneficial for the radial growth of larch in the Pangquangou area.